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We investigate a nonequilibrium coarsening dynamics of a one-dimensional Ising spin system with 
chirality. Only spins at domain boundaries are updated so that the model undergoes a coarsening 
to either of equivalent absorbing states with all spins + or — . Chirality is imposed by assigning 
different transition rates to events at down (H — ) kinks and up ( — h) kinks. The coarsening is 
characterized by power-law scalings of the kink density p ~ t^^ and the characteristic length scale 
^ ~ t^^^ with time t. Surprisingly the scaling exponents vary continuously with model parameters, 
which is not the case for systems without chirality. These results are obtained from extensive Monte 
Carlo simulations and spectral analyses of the time evolution operator. Our study uncovers the 
novel universality class of the coarsening dynamics with chirality. 

PACS numbers: 02.50.Ey, 05.50.+q, 05.70.Ln 



Coarsening takes place in various systems such as mag- 
netic systems, binary alloys, and social systems with 
opinion dynamics. When a system is quenched from a 
high-temperature disordered phase to a low-temperature 
ordered phase, a typical size of domains grows in time 
following a power law 



(1) 



with dynamic exponent z. It is known that coarsening 
systems are classified into a few universality classes de- 
pending on spatial dimensionality, order parameter sym- 
metry, conservation in dynamics, and so on 

The Ising model is one of the best studied coarsen- 
ing systems. It is symmetric under the global spin in- 
version {Z2 symmetry) and has a scalar order param- 
eter. Under the singlc-spin-flip Glauber dynamics Q 
that does not conserve the order parameter, the coars- 
ening dynamics is characterized by z = 2. On the other 
hand, the dynamic exponent is given by z = 3 under the 
spin-exchange Kawasaki dynamics Q conserving the or- 
der parameter. Systems with nonscalar order parameter 
constitute distinct universality classes [J. 

A coarsening process is rather simple in systems with 
discrete symmetry and nonconserving dynamics in one 
dimension. Consider a one-dimensional (ID) Ising spin 
chain with the Glauber dynamics at zero temperature, 
or equivalently the voter model In this model, only 
spins at domain boundaries can flip so that domain walls 
diffuse and annihilate in pairs. The diffusive nature sug- 
gests that the dynamic exponent is given by z = 2 and 
that the domain wall density decays algebraically as 



-.5 



(2) 



with an exponent S = 1/2. These scaling laws are verified 
by the exact solution [2, 0, 0| • 

The power-law scaling with 2 = 2 and 6=1/2 seems 
to be robust in one dimension. The (/-state Potts model 



with the zero-temperature Glauber dynamics exhibits the 
same scaling behavior 0, 0-1^. It is also observed in 
nonequilibrium systems. Consider the voter model with 
an additional exchange process of neighboring spins [lo| . 
It is equivalent to the branching annihilating random 
walk (BAW) model where domain walls diffuse, 

annihilate in pairs, and branch two offsprings. Despite 
the branching, the model displays the coarsening with 
the same exponents 12 , l3| ■ The voter model with a ki- 



netic constraint also displays the same scaling behavior 
with a logarithmic correction |14 1 . 



Most studies on the coarsening have focused on the 
role of order parameter symmetry On the other 

hand, some dynamical systems arc characterized by cou- 
pled symmetry and little is known about the coarsening 
dynamics in such systems. In this Letter, we investi- 
gate the coarsening dynamics of a ID Ising spin system 
which is invariant under the simultaneous inversion of 
spin and space. Remarkably, the model with the coupled 
symmetry constitutes a novel universality class that is 
characterized by continuously varying exponents. 

Our study was initially motivated by a flocking phe- 
nomenon of self-propelled particles. Flocking means a 
collective motion of particles, which is easily observed 
in a group of birds, fish, insects, or animals in nature. 
Vicsek et al |15| proposed a simple model where the mo- 
tion of each particle i at position is described by a 
velocity vector Vi of a constant speed. Each time step, 
the direction of Vi is updated to the average direction of 
particles within a fixed distance perturbed by a random 
noise. The system coarsens into a flocking phase when 
the noise strength is small 15, Note that spatial 

isotropy is broken spontaneously due to the motion of 
particles. Consequently, the model is symmetric under 
the simultaneous rotation/inversion of the velocity and 
the space. 

In order to investigate the effect of coupled symmetry. 
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FIG. 1. Plots of effective exponents 5eS for (a) u = 0, (b) 0.3, 
(c) 0.6, and (d) 0.9. 



we study an Ising spin system {cr = (ui, • • • , ajv)} with 
iji = ±1 in a ID lattice of N sites under periodic bound- 
ary conditions. Only spins near domain walls or kinks 
are updated following the rule 



- + 



u/2 



v/2 
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where + (— ) stands for a spin with value +1 (—1) and 
parameters over the arrows denote the transition rates 
of corresponding events. Without losing generality, we 
set u + V = 1 and u + v < 1. The model may be repre- 
sented with two kinds of domain walls or kinks: A and B 
corresponding to a down kink H — and an up kink — (-, re- 
spectively. The spin dynamics are then translated to an 
unbiased hopping of A (B) with rate v {v) and branch- 
ing of A ABA {B -)■ BAB) with rate u (u). Opposite 
kinks annihilate in pair upon collision. 

Our model is characterized by chirality: The transi- 
tion rates for events associated with down and up kinks 
are different. This chirality breaks the Z2 symmetry, 
but leaves system invariant under the simultaneous in- 
version of spin and space, cTj; — > — cr_i. Emphasizing 
the role of the chirality, the model will be referred to 
as the nonequilibrium chiral Ising model (NCIM). We 
remark that chirality is irrelevant for equilibrium Ising 
However, 



systems 
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it turns out to result in an inter- 
esting feature in nonequilibrium cases. 

The NCIM reduces to the voter model when v = v 
and u = u = 0, and the asymmetric simple exclusion 
process (ASEP) Q when v = v = 0. A mixture of them 
was studied in Refs. [H, [ill and was found to display 
complicated scaling behaviors. When u = u and v = v 
{Z2 symmetric case without chirality), discerning A from 
B becomes nominal and the model becomes equivalent to 



the BAW model 12|. It is solvable exactly [l3| and the 
kink density decays with (5 = 1/2 for all u < 1. Note 



TABLE I. Numerical values of 5 and z of the MCM for various 
values of it. For z, we present the results from the eigenspec- 
trum analysis [z (spectrum)] and from the seed simulations 
[z (seed)]. The numbers in parentheses indicate errors of the 
last digits. 



u 


5 


z (spectrum) 


z (seed) 





0.5000(1) 


2.0000(6) 


1.998(8) 


0.1 


0.4678(2) 


1.8789(5) 


1.877(9) 


0.2 


0.4346(3) 


1.7683(6) 


1.769(4) 


0.3 


0.4001(2) 


1.6666(6) 


1.667(2) 


0.4 


0.3639(5) 


1.5716(8) 


1.570(6) 


0.5 


0.3254(5) 


1.4818(7) 


1.483(4) 


0.6 


0.2837(4) 


1.3958(4) 


1.397(4) 


0.7 


0.2376(8) 


1.3117(8) 


1.313(6) 


0.8 


0.1850(8) 


1.227(2) 


1.228(3) 


0.9 


0.1195(3) 


1.136(5) 


1.139(6) 



that the NCIM is different from the so-called directed 
Ising model in which kinks are biased to a preferred di- 
rection 

First, we study the maximum chirality case with u = 
■D = where up kinks {B) are immobile and do not branch 
offsprings. This case is analogous to a ID version of the 
flocking model [isj if + (— ) spin is interpreted as a bird 
moving to the right (left) and if a bird is assumed to 
interact with another bird along its moving direction. 
For convenience, we refer to this case as the maximum 
chirality model (MCM). 

The system is prepared in an anti-ferromagnetic state 

(•••H i — •••) initially. Then, we measure the total 

kink density and average over Ns samples to obtain pit) 
for t < 10^. Just like the BAW model, we expect that 
p{t) decays in a power-law fashion for all w < 1 [2^ . 
So we investigate the behavior of an effective exponent 
defined as 



S,B{t) = -l0g[p{t)/p{t/b)]/l0gb 



(4) 



with a constant b. If the long time behavior of the kink 
density is given by p{t) = t^^{a + ct^'') up to a leading 
correction to scaling, the effective exponent in the long 
time limit behaves as Scs{t) o± S + ait~'^ , where oi is a 
constant depending on b. So when we draw 5eS against 
with the correct value of (, Scs{t) should approach 
to S with a finite slope as t — > 00. 

Figure [T] presents the behavior of the effective expo- 
nents for M = 0, 0.3, 0.6, and 0.9. The number of runs 
for all u is Ns = 5000 and system sizes for u = 0, 0.3, 
0.6, and 0.9 are N = 2'^*, 2^3, 2^^, and 2^1 , respectively. 
These system sizes are large enough that all runs did not 
fall into absorbing states where all spins have the same 
sign. 

The correction-to-scaling exponent ^ are roughly esti- 
mated from a fitting of S^s to the from S + ait~^ with 
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FIG. 2. Plots of 5 vs u (symbols) and its fitting function 
a2{l — u)^ with a2 ~ 0.5 and x ~ 0.62 (lines). Inset: Log-log 
plot of 5 against 1 — u together with the fitting function. 



three fitting parameters S, ai, and (. This procedure can 
be error-prone and the accuracy of the estimated C may 
be questionable. Still, however, if we draw i5eff against 
with C from a fitting, we observed that Jeff's intersect the 
y-axis with finite slope, which allows for estimating S. We 
observed that the corrections to scaling become stronger 
as u approaches 1, which is thought to be a signature of 
the crossover from the mean-field voter dynamics to the 
MCM; see below. The analysis of the effective exponent 
for other values of u is presented in Supplemental Mate- 
rial (2^ . Table |I] summarizes numerical results of 5 for 
various m's and Fig. [5] illustrates S against u. 

The intriguing feature of the MCM is that the expo- 
nent S varies with u. Furthermore, as u approaches 1, S 
seems to show a non-trivial power-law behavior. Indeed, 
if we fit 6 for the region u > 0.6 using 5 « 02(1 — u)^, 
with two fitting parameters 02 and x, we found that it 
fits the data quite well with 02 ~ 0.5 and x ~ 0.62; see 
Fig. [21 This singular behavior of 5 at w = f can be un- 
derstood as a result of a crossover from the mean-field 
voter dynamics to the MCM. 

When u is very close to 1 but not exactly I, the voter 
dynamics which occurs with rate v = 1 — u happens after 
many attempts of the ASEP dynamics. Since the sta- 
tionary state of the ASEP is totally uncorrelated Q , the 
voter dynamics can happen only after all spins are dis- 
tributed almost randomly. Hence, if we rescale the time 
as r = (I — u)t and take a limit u — > I with r kept fi- 
nite, the MCM should be the same as a mean-field voter 
model on a complete graph. Since coarsening does not 
occur on a complete graph of infinite size, S should be 
zero in the above mentioned limit. Thus, there should 
be a crossover from the mean-field voter dynamics to the 
ID MCM at w = 1. Note that the above argument does 
not depend on whether v and/or v are zero or not once 
V approaches zero as w — > I. Although the accuracy of x 



may be questionable, we can still insist that the crossover 
is described by the (non-trivial) exponent x- 

We substantiate the Monte Carlo results by study- 
ing the spectrum of the time evolution operator of the 
MCM. In general, a master equation can be mapped to an 
imaginary time Schrodinger equation with Hamiltonian 
H whose eigenvalues contain most of relevant information 
of the system. For instance, the directed percolation sys- 
tem has been studied successfully with the eigenspectrum 
analysis 23l425|. 

For the MCM, the Hamiltonian takes the form 

^ N N 

H = T ^(1 + -^.f )(1 ' ^-+1) -uY,K<^t+i 



(l-u) 



N 



^{a-(l-af+i) + (l+<7,^)<7+ J,(5) 



where ai is the Pauli spin operator acting on a spin at site 
i. We label eigenvalues of H as En with n = f , • • • ,2^ 
and call En the energy of the n-th level. Since H is not 
Hermitian, En may have a complex value. The eigenval- 
ues are sorted in the ascending order of 3fi[i5„], the real 
part of En ■ There are two trivial levels with Ei = E2 ^ 
corresponding to the two absorbing states with all spins 
having the same sign. Other low-lying energy levels with 
n > 2 define the relaxation time as r„ — 1/(3? [£'„]). 

We have diagonalized numerically the Hamiltonian up 
to iV = 20 to obtain the longest relaxation time T3 . Since 
T3 ^ , the dynamic exponent z is estimated by extrap- 
olating an effective exponent Zoff(iV) = \n[T3{N) / T3{N — 
2)]/ln[N/{N - 2)] with the Bulirsch-Stoer (EST) algo- 
rithm 0, [i^l . The results are summarized in Table HI 
As an illustration, detailed analysis of Zeff for u = 0.8 
along with the extrapolation according to the EST algo- 
rithm is presented in Supplemental Material [22 1. 

We have also performed independent Monte Carlo sim- 
ulations starting with a single down kink A in the middle 
of an infinite lattice, which are generally referred to as 
seed simulations. We have measured the particle spread- 
ing distance ^{t) ~ t^^^ to obtain the dynamic exponent 
z following the similar effective exponent analysis as done 
for S. Both results for z are in perfect agreement with 
each other and vary continuously with u. See Supple- 
mental Material [22| for the full analyses of z. 

We have shown that the coarsening dynamics of the 
MCM is characterized by continuously varying critical 
exponents. Note that chirality lies both in the ASEP 
events (u ^ u) and the voter- model events {v 7^ v). In 
order to investigate which one is the essential ingredient, 
we studied two more cases: One is the case with u = u 
and V V = which will be called the symmetric exclu- 
sion and chiral voter (SECV) model and the other is the 
case with u ^ u ^ and v = v to he called the chiral ex- 
clusion and symmetric voter (CESV) model. The CESV 
model is a particular limiting case of the model studied 
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1/lnt 

FIG. 3. (Color online) Plots of S^fi as functions of 1/lnt 
for u = 0.1, 0.2, 0.5, and 0.8 (from top to bottom) for the 
SECV. A line segment with slope —0.5 is for a guide to the 
eyes. Inset: Plot of p{t)Vt/{lnt)°-^^ vs t for the SECV with 
M = 0.8 on a semi-logarithmic scale. 



in Ref. 0, [l^. Remind that v is always set to I — u. 

The analyses otScS for the SECV with u = 0.1,0.2,0.5, 
and 0.8 are summarized in Fig. [31 It seems that 6cs{t) for 
all u approaches 1/2 with logarithmic corrections. Notice 
that if there exists a logarithmic correction as 



p{t) - (lnt + C)«A*' 



(6) 



with a constant C, the effective exponent defined in 
Eq. (HI) should behave as 6ee{t) sa 6 — K/\nt in the asymp- 
totic regime. Thus, if we plot Scff{t) as a function of 
1/ln^, the effective exponent should intersect the y-axis 
with slope —K. This phenomenon is quite pronounced 
for the cases oi u = 0.5 and 0.8 and the slope seems 
to be around 0.5. Indeed, if p{t)\/t/{\nt)°-^^ is plotted 
against t on a semi-logarithmic scale (see Inset of Fig. [3] 
for the case of it = 0.8), a flat region is observable in 
the long time limit for more than two log-decades. Al- 
though the accurate value of n is hard to estimate, we 
can conclude that there exist a systematic logarithmic 
correction as shown in Eq. ^ with the leading scaling 
exponent S — 1/2 unchanged in the SECV model. 

A logarithmic correction in a ID coarsening has been 
reported in a different model |14| that corresponds to 
the voter model with a weak kinetic constraint. In 
that model, the kink density decays faster than as 
l/(Vtlnt), but in our case it decays slower than 1/Vt. 
Qualitatively, this slowing down should be attributed to 
the presence of branching dynamics which increases the 
number of kinks. However, a quantitative analysis re- 
quires further investigation, which is beyond the scope 
of this Letter. For our purpose, it is enough to conclude 
that the continuously varying exponents are not due to 
the chiral voter dynamics. 





FIG. 4. Plots of (5cff as functions of l/ln(t) for u = (a) 0.2, 
(b) 0.5, (c) 0.6, and (d) 0.8 for the CESV. For the case of it = 
0.5 and 0.6 [(b) and (c)], fitting results of a function 5 + d/lnt 
with 5 and d fitting parameters are also depicted. 



The CESV model shows a more intriguing feature. 
We present the effective exponent data for the density 
decay in Fig. 31 For u = 0.2 [Fig. HJa)], 6cs seems 
to approach 0.5 with negligible logarithmic correction. 
For u = 0.5 [Fig. HJb)], we cannot make a firm con- 
clusion whether S < 0.5 or (5 = 0.5 due to strong 
correction-to-scaling behavior. Quite interestingly, when 
u > 0.5 [Fig. lU^c) and (d)], 6 deviates from 0.5 signifi- 
cantly even under the assumption of a logarithmic cor- 
rection. So we conclude that the CESV has continuously 
varying exponents with possible logarithmic corrections 
when u > 0.5. This study shows that the chirality in the 
spin exchange is responsible for the continuously varying 
exponents. Nevertheless, it remains open why and when 
there appears a logarithmic correction in the coarsening 
process. 

To summarize, we have studied the one-dimensional 
coarsening dynamics of nonequilibrium Ising spin sys- 
tems with chirality. Although chirality is irrelevant in 
equilibrium Ising systems, it turns out that the chirality 
can lead to continuously varying scaling exponents in the 
nonequilibrium chiral Ising model. In particular, it turns 
out that the chirality in spin exchange plays a crucial 
role. 

It is rare to observe continuously varying exponents 
from systems without quenched disorder. The g-state 
Potts model with zero-temperature Glauber dynamics 
was studied in Refs. [1, It was found that the crit- 
ical exponent describing the power-law decay of the per- 
sistent probability varies continuously with q. Neverthe- 
less, the coarsening dynamics is still pure diffusive and 
characterized by z = 2 and (5 = 1/2 at all values of q. 
We notice that continuously varying exponents were re- 
ported in a ID sandpile model without dissipation (28j 
and that there is actually a parallelism between this sand- 
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pile model and the MCM. This connection will be dis- 
cussed elsewhere 23. 



Some of Isinff spin systems are exactly solvable in 
one dimension [J [23| , for equations governing the 
time evolution of correlation functions are closed. In the 
presence of the chirality, however, the equations are not 
closed, which makes the exact solution for the NCIM not 
available in general. Our numerical finding of the uni- 
versality class with continuously varying exponents can 
be established more firmly if one find a minimal con- 
tinuum equation obeying the proper symmetry property. 
We leave it as a future work [29[. 
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